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In this work we derive the modifications to the Boltzmann equation governing the cosmic evo- 
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^ ■ I. INTRODUCTION 

m 

In a previous work [1| we have discussed a case study of dissipative Liouvillc-string cosmology, involving non 
critical string cosmological backgrounds, with the identification Q of target time with the world-sheet zero mode of 
the Liouville field [3j- Such cosmologies were found to asymptote (in cosmic time) with the conformal backgrounds 
of 3, which are thus viewed as equilibrium (relaxation) configurations of the non-equilibrium cosmologies. 

In terms of microscopic considerations, a model for such departure from equilibrium could be considered the collision 
ij of two brane worlds, which results [SJ in departure from conformal invariance of the effective string theory on both 
the bulk and the brane world, and thus in need for Liouville dressing to restore this symmetry |3|]. Dynamical 
arguments 0, then, stemming from minimization of effective potentials in the low-energy string-inspired effective 
field theory on the brane, imply the (eventual) identification of the zero mode of the Liouville mode with (a function) 
of target time There is an inherent time irreversibility in the process, which is associated with basic properties of 
the Liouville mode, viewed as a local (dynamical) renormalization group (RG) scale on the world-sheet of the string. 
This implies relaxation of the associated dark energy of such cosmologies, and a gravitational friction, associated with 
the conformal theory central charge deficit. For a recent review we refer the reader to @, where concepts and methods 
arc outlined in some detail. 

For our purposes in this letter we would like to concentrate on one interesting aspect of the non-critical string 
cosmologies, associated with the off-equilibrium effects on the Boltzmann equation describing relic abundances and 
the associated particle-physics phenomenology. Indeed, in conventional cosmologies, a study of relic abundances by 
means of Boltzmann equation that governs their cosmic time evolution yields important phenomenological constraints 
on the parameters of particle physics (supersymmetric) models using recent (WMAP |7j and other) astrophysical 
Cosmic Microwave Background (CMB) data. Essentially, the astrophysical data constraint severely in some cases 
the available phase-space distributions of favorite supersymmetric dark matter candidates such as neutralinos 0. 
When off-equilibrium, non-critical string cosmologies are considered []] , which notably are consistent with the current 
astrophysical data from supernovae, as demonstrated recently 0, then there are significant modifications to the 
Boltzmann equation, stemming from extra sources (dilaton) and off-shell (non-critical, non equilibrium) terms, which 
affect the time evolution of the phase-space density of the species under consideration. 

It is the purpose of this paper to derive such modifications in detail, and then use them in order to discuss briefly 
particle-physics models constraints, especially from the point of view of supersymmetry. With regards to this last 
issue, in this article we shall present an approximate analytical treatment, which will only provide hints to what may 
actually happen. A complete analysis requires numerical studies which are postponed for a future publication. 
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II. MODIFIED BOLTZMANN EQUATION IN NON-CRITICAL STRINGS 

Consider the phase space density of a species X, which is assumed coupled to the off-shell (non-critical string) 
background terms: 

f(\p\,t) 

where quantities refer to the Einstein frame, and the Einstein metric is assumed to be of Robertson- Walker (RW) 
type. Throughout this work we follow the normalization and conventions of 

For completeness we state the main relationships between the Einstein and cr-model frames Q: 

9v» = e S M „ , ^ = e (2.1) 

where $ is the dilaton field, and the superscript a denotes quantities evaluated in the a-model frame. 

Since in the RW Einstein frame goo = — 1: and in the cosmic Einstein co-moving frame we have for a generic massive 
species with mass m (that we shall be interested in in this work) p£ = mdx^/dr = m(dx^ / dt)(dt / dt a ) = mp^e - *, it 
can be readily seen that 

/ 3 \ V 2 

\p a \ = (I''-'' I^UJX 1 = (pV.%)' /2 = \P\ = a(t) ( 5>¥j (2.2) 

We shall be interested in the action of the relativistic Liouville operator L, acting on a phase-space density of 
species X, f(\p\,t,g l ), which in general depends on the off-shell backgrounds g % = (gu,$), where gu are the spatial 
components of the RW space-time metric in Einstein frame, and $ is the dilaton. The off-shell dependence comes 
about due to the interpretation of target time as a local world-sheet renormalization-group scale, the Liouville field. 

We commence our analysis by recalling the relativistic form of the Liouville operator in conventional general 
relativity: 

The second term on the r.h.s. is the relativistic form of the force, following from the geodesic equation. For a RW 
Universe, only the time-energy part survives from the first term, i.e. 

*"»t-*5- (2 ' 4) 

Moreover, the connection part receives non trivial contributions only from the terms J] i T^p l p 1 -^, i = 1,2,3a spatial 
index (assuming for concreteness an already compactified string theory, or a theory on a three-brane world), with 
r° = —ad, where the overdot denotes derivative w.r.t. cosmic RW time t, identified in our string theory with the 
Einstein- frame RW time (|2.1|) . 

Hence, the conventional part of the Liouville operator in a RW Universe would read [Ifl |: 

v - E m - aa T,py^ = E oi--X ^im (2 - 5) 

i— 1 i— 1 

The presence of off-shell, non-critical (Liouville) string Q backgrounds, upon the (dynamical) identification Q,[ME1 
of the Liouville mode with the target time, one will receive extra contributions in the expression for the associated 
Liouville operator, stemming from the fact that the latter is nothing but a total time derivative. 

In the context of our string discussion it is important to specify that the initial frame, from which we commence 
our discussion, is the so-called Einstein frame. It is in this frame that the usual RW cosmology is obtained in string 
theory Q, for which the expression l|2.5|l is valid. 

To discuss the non-critical-string (off-shell) corrections to Boltzmann equation, it will be necessary to consider time 
derivatives in the a-model frame. This is due to the fact that it is in this frame that the target time X° = t a is 
related simply to the Liouville mode tp in non-critical string theories 0, 



<y3 + t a = 



(2.6) 
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This relation is obtained dynamically from minimization arguments of the low-energy effective potential of some 
physically interesting cosmological models, for instance those involving colliding brane worlds 0,0- To be precise, 
the initial relation l|2.6|) derived in the specific model of |(| reads: ip/y/2 + t a = , with the factor of y/2 arising 
from (logarithmic) conformal-ficld-theory considerations, and being crucial for yielding a Minkowski space-time in 
the scenario of 6], where the coordinate X° assumed an initial Euclidean signature, appropriate for a world-sheet 
path-integral quantization. Nevertheless, in our approach below we shall absorb the \[2 into the normalization of the 
units of the Regge slope a' , for convenience. 

In such scenarios the target-space dimensionality of the string is extended to D + 1 initially, with two time-like 
coordinates, t and (p. Eventually, these two coordinates are identified (c.f. I|2.6[l ) in the solution of the generalized 
conformal invariance conditions that express the restoration of the conformal invariance by the Liouville mode, as 
we discuss below. This implies that, initially, the Liouville operator l|2.3[l acquire extra Liouville components (^p* 3 ), 
with p v = E v . Hence, there are additional structures on the right hand side of (|2.3fl . of the form E v -^-, together with 
the corresponding amendments in the Christoffcl-symbol dependent d/dp^ terms, stemming from the extra Liouville 
coordinate on the target space time, which implies appropriate extra components of the extended-target-space-time 
metric. The (dynamical) implementation of l|2.6|l restricts oneself to a /^-dimensional hypersurface in the extended 
target-space, with the identification of E v — > E, the ordinary (physical) energy. From now on we restrict ourselves 
on this hypersurface, bearing in mind however that the identification of the Liouville mode tp with the target time 
should only be implemented at the end, and thus any dependence on (p should be kept explicit at intermediate steps 
of the pertinent calculations. 

The connection between time derivatives in the Einstein and stringy-cr-model frames is provided by the chain rule 
of differentiation = using l|2.1|l . This will result in multiplicative factors of e~* in front of the appropriate 

non-critical-string modifications of the Liouville operator l|2.5[) . 

Since, as mentioned above, the Liouville operator is essentially a total time derivative operator, and in our case time 
is related to a world-sheet renormalization group scale, p(cr, r), taken to be local on the two-dimensional surface as a 
result of world-sheet general covariance p"H , the sought-for non-critical-string modifications of the Boltzmann equation 
emerge from the implicit dependence of the g l background fields on p(cr, t), that is the corresponding /3-functions, 
which in a critical-string theory would vanish. 

We now remark that in the approach of 0], the local RG scale /?(£) was identified with the dynamical Liouville mode 
tp, which implied that in this approach the local conformal invariance of the world-sheet of the string was restored 

p = p. (2.7) 

With these in mind we then modify the relativistic form 12. 3[) . (|2.4|) by replacing 

d d _ d f dt a d P (0 dg' H " d 



dt dt dt 7 S dt dt a d P (0 dg' H " 

where £ = cr, r, denotes the world-sheet coordinates, J s is a world-sheet integration, and the index u i" runs on both, 
a (discrete) background field space, {gu(y), §{y)}, i = 1,2,3 and a continuous D-dimensional (target) space-time y. 
Hence, summation over "i" includes integration J d D y^/—g(...) . Note that since the (physical) target time in the 
cr-model frame, t a , is only the world-sheet zero mode of the corresponding cr-model field X°(£), eventually the world- 
sheet integration J" E in (|2.8|1 disappears, since only the world-sheet zero mode p of the local RG scale p(£) (Liouville 
mode) will yield a non zero contribution. In the last equality of the right-hand side of 12.8fl . the quantity r\ is defined 
as: r\ = Jjp-. If we use the identification p(£) = cp(£) and Eq. H2.6|l . then rj = —1, but at this stage we keep it general 
to demonstrate explicitly the renormalization-scheme dependence (i.e. choice of the local scale p(£))- 

The non-critical string contributions ^p- = f3 z are the Weyl anomaly coefficients of the string, but with the renor- 

malized cr-model couplings g l being replaced by the Liouville-dressed [3j quantities. This is a particular feature of 
the approach of 0] viewing the Liouville mode as a local world-sheet RG scale, as mentioned previously. In turn the 
latter is also identified (c.f. (|2.6|) ') with (an appropriate function of) the target time t. 

The dynamics of this latter identification is encoded in the solution of the generalized conformal invariance condi- 
tions, after Liouville dressing, which read in the cr-model frame 

-p = /' + Q/ (2.9) 

where the prime denotes differentiation with respect to the Liouville zero mode p, and the overall minus sign on the 
left-hand side of the above equation pertains to supercritical strings Q, with a time like signature of the Liouville 
mode, for which the central charge deficit Q 2 > by convention. 
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Notice the dissipation, proportional to the (square root) of the central charge deficit Q, on the right-hand side 
of 12. 9[) . which heralds the adjective Dissipative to the associate non-critical-string- inspired Cosmological model. 
Moreover, the Weyl anomaly coefficients f3 l , i = {$, g^}, whose vanishing would guarantee local conformal invariance 
of the string-cosmology background, are associated with off-shell variations of a low-energy effective string-inspired 
action, 

= QaP Qij = Um z ^ Q z 2 z{Vi{z, z)Vj{0, 0)) (2.10) 

with z,z (complex) world-sheet coordinates, Qij the Zamolodchikov metric in theory space of strings, and V* the o~- 
model vertex operators associated with the a- model background field g l . It is this off-shell relation that characterizes 
the entire non-critical (Q) Cosmology framework, associated physically with a non-equilibrium situation as a result 
of an initial cosmically catastrophic event, at the beginning of the (irreversible) Liouville/cosmic-time flow |5(. 

The detailed dynamics of l|2.9|l are encoded in the solution for the scale factor a(t) and the dilaton <f> in the simplified 
model considered in pj, after the identification of the Liouville mode with the target time l|2.6|l . In fact, upon the 
inclusion of matter backgrounds, including dark matter species, the associated equations, after compactification to 
four target-space dimensions, read in the Einstein frame [l|: 

e 2 * 

3 H — Q m — = Qq, 

Qii 

2 H + g m + e$ + p m + p$ = — 

a z 

^ , o tt^ , 1 9V a u 1 _ — . 3 Qu e 2 * - . 

$ + 3H$ + - — + - (g m - 3p m ) = --—- —g 9 . (2.11) 

where g m (p m ) denotes the matter energy density (pressure), including dark matter contributions, and g$ (p$) the 
corresponding quantities for the dilaton dark-energy fluid. All derivatives in H2.11|l are with respect the Einstein time 
t which is related to the Robertson- Walker cosmic time taw by t — uj tuw Without loss of generality we have taken 
lj = \/SHq where Hq is the Hubble constant. With this choice for u> the densities appearing in (|2.11l) are given in 
units of the critical density. Then one can see for instance that if the time t^w is used the first of equations (|2.11|) 
above receives its familiar form in the RW geometry when the contributions of the dilaton and non-critical terms are 
absent. 

The overdots in the above equations denote derivatives with respect to the Einstein time. Their right-hand side 
contain the non-critical string off-shell terms: 

= e - 2 * (* - $ 2 + Q e *$) 
Gu = 2« 2 (f + 3ff$ + $ 2 + (1 - q)H 2 + Qe*(<i> + H) ) . (2.12) 

In the above equations H = (a) /a is the Hubble parameter and q is the deceleration parameter of the Universe 
q = —da/a 2 , and are both functions of (Einstein frame) cosmic time. The potential V a ii appearing above is defined 
by V a ii — 2Q 2 exp (2$) + V where, in order to cover more general cases, we have also allowed for a potential term 
in the four-dimensional action — J d 4 y^/—G V in addition to that dependent on the central charge deficit term Q. 
Although we have assumed a (spatially) flat Universe, the terms on the r.h.s., which manifest departure from the 
criticality, act in a sense like curvature terms as being non-zero at certain epochs. The dilaton energy density and 
pressure are given in this class of models by: 

= \ ( 2$ 2 + Van) 

Pi = ~ ( 2$ 2 - V a ii) ■ (2.13) 

Notice that the dilaton field is not canonically normalized in this convention and its dimension has been set to zero. 

For completeness, we mention at this point that the dependence of the central charge deficit Q(t) on the cosmic 
time stems from the running of the latter with the world-sheet RG scale 0, , and is provided by the Curci - Paffuti 
equation [lSj expressingthe renormalizability of the world-sheet theory. To leading order in an a' expansion, which 
we restrict ourselves in 1] and here, this equation in the Einstein frame reads: 



= _6e- 2 *(ff + $)^ 2 



6e - 2 * + (2.14) 
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For future use we also state here the corresponding continuity of the matter stress tensor, which is not an independent 
equation, but can be obtained from l|2.11[l by appropriate algebraic manipulations: 

^ + 3H(g m + p m ) + - $ (g m - 3p m ) = 6 (H + $) ■% • (2-15) 

This expresses the non-conservation equation of matter as a result of its coupling to both, the dilaton source terms 
and the off-shell, non-equilibrium (non-critical-string) backgrounds. 

A consistent solution of a(t), &(t), and the various densities, including back reaction of matter onto the space-time 
geometry, has been discussed in , where we refer the interested reader for further study. Also note that a preliminary 
comparison of such non-critical strings theories with astrophysical data, demonstrating consistency at present, is given 
in Q. 

After this necessary digression we now come back to discussing the derivation of the dilaton-source and non- 
critical-string induced modifications to the Boltzmann equation, governing the cosmic evolution of the various species 
densities. It is important for the reader to bear in mind already at this stage that the Boltzmann equation does 
not contain any independent information from the dynamical equations H2.II(I . but it should be rather viewed as 
an effective way of describing the cosmic evolution of the density of a given species, consistent with the continuity 
equation l|2.15|l for the total matter energy density. We shall come back to this important point later on. 

At the moment, let us concentrate first on the dilaton-source and non-critical-string background contributions to 
the Liouville operator l|2.8[) . The presence of (time-dependent) dilaton source terms implies an explicit dependence 
of the phase space density of a species / on <E>, while the non-conformal on the world-sheet) nature of the metric and 
dilaton background, induce a Liouville mode p dependence through the corresponding backgrounds: 

f(\0\,tMt,P),9v»&P)) ■ (2-16) 

The non-critical string terms can be expressed, as we have seen above, in terms of the corresponding Weyl anomaly 
coefficients, which are non zero as a result of departure from conformal invariance of the pertinent string background. 
Despite their non-linear looking appearance when expressed in terms of the Weyl anomaly coefficients (which depend 
on the Ricci tensor and second covariant derivatives of the Dilaton field), such terms acquire a particularly simple 
linear form once the identification of time with the Liouville mode (I2.6f) is implemented, which in effect implies that 
the explicit solution of l|2. 9JI . (|2.II(I must be taken into account when discussing the Boltzmann equation. 

With the above in mind, the form of the non-critical-string (off-shell) and dilaton-source corrections to the Liouville 
operator reads (expressing quantities in the Einstein frame, taking into account the fact that in the cr-model frame 
g" and <& are treated as independent field-theory variables, and using (|2.1|0 : 

L " U ' dil = E *d* + ET > e P W* 7 ^ E *d* + ^ (jg dgj dg- + d-pd$) = 
B*£ + ^ { *(± + «)±£ + *£ } (2.17) 

where the overdot denotes the derivative with respect to the (Einstein frame) cosmic time t, and the prime the 
derivative with respect to p, and we used the fact that: 

M = 2 a 2 e 2 * (- + $')» ^ = e~ 2 * (2.18) 
dp a dg^ 

with goo = — 1) gu = a 2 , an d ^ i = /3^ lav the corresponding Weyl anomaly coefficient, in order to express /3^ rav 
in terms of quantities in the Einstein frame. We have also been careful to separate the dilaton dependence on p(t) 
from that on the ordinary time t (prior to the identification p.fijl ). We remind the reader, at this point, that the 
identification 1)2. 6J) should only be applied at the end of the calculations, and should be viewed as a specific choice 
of "renormalization" scheme. The "renormalization" scheme l|2.6|l . i.e. rj = — 1, which is obtained by a physical 
principle (dynamics) in certain models of interest to us here may thus be viewed in this context as the "physical 
scheme" , where observational cosmology can be analyzed. In such a scheme there is a simplification of the complicated 
non-linear expressions encoded in the Weyl anomaly coefficients. 

/ 3 . . \l/2 

Taking into account the implicit dependence of /(|p|, t, $) on gu through \p\ = I ^2 i= iP l p l gu ) , it is evident that 
df/dgu = (d\p\/dgu) df/d\p\ = Jj4, from which we can easily see that the final form of the action of the Liouville 
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operator, including dilaton/source terms and non-critical string corrections, on the phase-space density / (12.16(1 . is 
such that: 



C' 



£off-shcll,dil) / — C[f] , 



§RfI-" e *(^ + *') - (* + -**') I + 1™ ■ <- j > 

Upon considering the action of the above operator on the density of a given species X, n = J d 3 pf we then arrive, 
after some straightforward momentum integration by parts, at the modified Boltzmann equation for a four-dimensional 
effective field theory after string compactification (or restriction on three-brane worlds), in the presence of non-critical 
(off-shell) string backgrounds and dilaton source terms: 

f - 1 / < 1 - (* + / *>& - (7 + *') / + / ^ - 

f + 3% - 3,e« £ + *) „ + / £cw - (4 + „•*) I ftg M) 

The collision term C[/] assumes the usual form in conventional particle cosmology [ljj. The reader is invited to 
compare the final equation for the cosmic time evolution of the density n in the second line of (|2.20() to the continuity 
equation for the total energy density of matter (|2.15|l . As discussed previously, the Boltzmann equation l|2.20|l should 
be compatible in the sense of leading to no extra information) with the conservation equation l|2.15|l , as well as the 
(modified) Einstein equations (|2.11|) . 

In the above scenario, where the non-critical string contributions have been obtained through the identification of 
time with the (world-sheet zero-mode of the) Liouville field p ((2.6(1 . there is some universality in the coupling of all 
matter species, including dark matter, to these off-shell background terms, which may be traced to the equivalence 
of all species coupled to gravity. This is to be contrasted with the more general phenomenological case studied in Q , 
where scenarios involving only the coupling of the non-critical string backgrounds with the exotic dark matter species 
have been considered as well. 

The final issue to be discussed in this section pertains to the form of the dependence of / on the dilaton source 
terms, which would survive a dilaton-driven critical-string cosmology case, such as the one considered in |13|. We 
constrain this form by requiring that in the Einstein frame there are two types of dependence on <£>: (i) explicit, of 
the form e~ 4 *, arising from the fact that in our approach, the phase space density is constructed as a quantity in 
the cr-model frame of the string, which is then expressed in terms of quantities in the Einstein frame. As such, it is 
by definition (as a density) inversely proportional to the proper cr-model frame volume V a — J d 4 Xy / —g' 7 oc e 4 * on 
account of ((2.1(1 ; (ii) implicit, corresponding to a dependence on $ through the Einstein- frame metric gu 1(2.1(1 . Hence 
the general structure of / is of the form: 



= e^V; i)cxe~ 4 *^(|pl,x,t) (2.21) 



This implies that: 



" d\p\ df 



-An- J d 3 p\p\^ = -4n + 3 J d'pf \\p\,t) = -n , (2.22) 

where in the last step we have performed appropriate partial (momentum-space) integrations. 
The final form of the Liouville operation 1(2.20(1 . then, reads: 



9$ dga J dga d\p\ 

;3 



, 3 - n - $n = Sr/e^—n + Arie m &n + / — C[f] (2.23) 
at \a J a J E 

We now notice that non-critical terms can be expressed in terms of the Weyl anomaly coefficients for the (cr-model) 
graviton and dilaton backgrounds as: 

3^=4^-* (lg^ff r _p$ e 2* 

4^$' = 477e* /3* (2.24) 
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where we used the Einstein frame metric to contract indices, with /3 Grav denoting the graviton Weyl anomaly coeffi- 
cient. In H'2.24f) we have taken into account that in the class of models of 0,^3, we are concentrating in this work, 

/3 G o av = 0. 

Thus, we observe from l|2.23[l . that terms proportional to the dilaton /?* cancel out, leaving only the graviton 
contributions as the non-critical string corrections to the Boltzmann equation, which finally becomes: 

d i + 3 (:) " - *» - b (» - v-#r) » + / %w (2.25) 

It is important to recall once more that, upon Liouville dressing, which restores the conformal invariance of the model, 
the graviton and dilaton Weyl anomaly coefficients satisfy 1)2. 9fl . which upon the identification of the Liouville mode 
with (a function of the cosmic) time yield p.IIfl • 

We next remark that in the case with non-critical string backgrounds, consistency requirements between the two 
equations determine the back reaction effects of matter ont o sp ace time. For instance, in the absence of matter, in the 
specific (but generic enough) non-critical string model of [l4|, we have that the present-era dilaton and scale factor 
of the Universe are such that Q ~ 0, $ + 1=0. In this particular case the off-shell terms at late eras are practically 
zero, and the system has reached equilibrium. The corresponding Boltzmann equation for the dilaton species densities 
is then satisfied in agreement with the corresponding continuity equation l|2.15|) in absence of matter. 

However, in the presence of matter, we know from the analysis of Q that 4> + - ^ 0, and the right-hand-side of 
the Boltzmann equation l|2.20|l . with n now representing matter species, yields terms proportional to n. Consistency 
with the continuity equation l|2.15|l , then, implies that the off-shell (non-critical string background) terms in the latter 
should be proportional to the total number density of matter species (back reaction effects). This and other similar 
consistency checks should be made when considering solving such non-equilibrium cosmologies. 



III. PHENOMENOLOGY OF PARTICLE PHYSICS MODELS AND MODIFIED BOLTZMANN 

EQUATION 

In this section we consider solutions of the modified Boltzmann equation (|2.23|) , or equivalently (|2.25|l , for a particle 
species density n in the physically interesting case of supersymmetric dark matter species, such as neutralinos, viewed 
as the lightest supersymmetric particles (LSP). Such Cold Dark Matter candidates lead to a rich phenomenology of 
supersymmetric particle physics models. In the context of conventional Cosmology Q, some of these models can be 
constrained significantly by the recently available astrophysical data on Cosmic Microwave Background temperature 
fluctuations p|. The calculation of relic abundances will be done in some detail, in order for the reader to appreciate 
better the role of the non conventional terms in (|2.25() . . 

It is convenient to write the Boltzmann equation for the density of species n in a compact form that represents 
collectively the dilaton-dissipative-source and non-critical-string contributions as external-source T(t)n terms: 

f + = T(t)n + J ^C[f] , T(t) ee $ + \ne~*^^ (3.1) 

where we work in the physical scheme <|2.6ll from now on, for which r\ = — 1. Depending on the sign of T(t) one has 
different effects on the relic abundance of the species X with density n, which we now proceed to analyze. To find an 
explicit expression for T(t) in our case we should substitute the solution of l|2.9|) . more specifically (|2.11|l . analyzed in 
. Regarding the form of i|3.I[l it is nice to see that the extra terms can be cast in a simple-looking form of a source 
term T(t)n including both the dilaton dissipation and the non-critical-string terms. Of course T{t) is complicated and 
requires the full solution outlined in 

In a more familiar form, the interaction term C[f] of the above modified Boltzmann equation can be expressed in 
terms of the thermal average of the cross section a times the Moeller velocity v of the annihilated particles 01 



i: = -3-n-(va) (n 2 - n%) + T n (3.2) 
at a H 

Before the decoupling time tf, t < tf, equilibrium is maintained and thus n = n eq for such an era. However, it is 
crucial to observe that, as a result of the presence of the source T terms, n eq no longer scales with the inverse of the 
cubic power of the expansion radius a, which was the case in conventional (on-shell) cosmological models . 
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To understand this, let us assume that n = n^q at a very early epoch t$. Then the solution of the modified 
Boltzmann equations at all times t < tf is given by 

n eq a 3 = n$ a 3 {t ) exp ( / Tdt) . (3.3) 



The time to characterizes a very early time, which is not unreasonable to assume that it signals the exit from the 
inflationary period. Soon after the exit from inflation, all particles are in thermal equilibrium, for all times t < tf, 
with the source term modifying the usual Boltzmann distributions in the way indicated in Eq. (|3.3I) above. It has 
been tacitly assumed that the entropy is conserved despite the presence of the source and the non-critical-string 
contributions. In our approach this is an approximation, since we know that non-critical strings lead to entropy 
production. However, as argued in our previous works on the subject 0, the entropy increase is most significant 
during the inflationary era, and hence it is not inconsistent to assume that, for all practical purposes, sufficient for our 
phenomenological analysis in this work, there is no significant entropy production after the exit from inflation. This 
is a necessary ingredient for our approach, since without such an assumption no predictions can be made, even in the 
conventional cosmological scenarios. Thus, the picture we envisage is that at to the Universe entered an equilibrium 
phase, the entropy is conserved to a good approximation, and hence all particle species find themselves in thermal 
equilibrium, despite the presence of the T source, which slowly pumps in or sucks out energy, without, however, 
disturbing the particles' thermal equilibrium. 

From the above discussion it becomes evident that it is of paramount importance to know the behaviour of the 
source term at all times, in order to extract information for the relic abundances, especially those concerning Dark 
Matter, and how these are modified from those of the standard Cosmology. Before embarking on such an enterprize 
and study the phenomenological consequences of particular models predicting the existence of Dark Matter, especially 
Supersymmetry-based ones, we must first proceed in a general way to set up the stage and discuss how the relic density 
is affected by the presence of the non conventional source terms discussed above. 

For the sake of brevity, we shall not deploy all the details of the derivation of the relic density, but instead 
demonstrate the most important features and results of our approach, payi ng p articular attention to exhibiting the 
differences from the conventional case. Generalizing the standard techniques jlOj . we assume that above the freeze-out 
point the density is the equilibrium density as provided by Eq. I|3.3|l . while below this the interaction terms starts 
becoming unimportant. It is customary to define x = T/m% and restrict the discussion on a particular species x OI 
mass m^, which eventually may play the role of the dominant Dark Matter candidate. It also proves convenient to 
trade the number density n for the quantity Y = n/s, that is the number per entropy density 10] . The equation for 
Y is derived from l|3.2(l and is given by 

- = m, (va) ( -G N g eff ) (*+-_) (Y - Y? q) - - (1 + --) Y . (3.4) 

where Gn = ^-/Mp lanch is the four-dimensional gravitational constant, the quantity H is the Hubble expansion rate, 
h denote the entropy degrees of freedom, and (va) is the thermal average of the relative velocity times the annihilation 
cross section and g e f t is simply defined by the relation |ic| 

g + Ag = ^T 4 ~g eff . (3.5) 

The reader should notice at this point that Ag incorporates the effects of the additional contributions due to the 
non-critical (off-shell) terms and the dilaton dissipative source, which are not accounted for in the g e // of conventional 
Cosmology |lfl |. hence the notation g e f f ■ We next remark that p, as well as Ap, as functions of time are known, once 
one solves the cosmological equations. However, only the degrees of freedom involved in p are thermal, the rest, like 
the comological-constant term if present in a model, are included in Ap. Therefore, the relation between temperature 
and time is provided by 

P=^ T4 9eff(T) (3.6) 

while p + Ap are involved in the evolution through (c.f. I|2.11|l ) 

H*= 8 -^L(p + A P ) . (3.7) 

Thus, it is important for the reader to bear in mind that Ap contributes to the dynamical expansion, through Eq. 
(13.71) . but not to the thermal evolution of the Universe. The quantity g e ff, defined in (|3.5|l . is therefore given by 

30 

9eff = 9eff + ^T^Ap . (3.8) 
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The meaning of the above expression is that time has been replaced by temperature, through Eq. I|3.6[l . after solving 
the dynamical equations. In terms of g e ff the expansion rate H is written as 



H 2 = 



47r 3 Gjv 4 . 

-45- *// 



(3.9) 



This is used in the Boltzmann equation for Y and the conversion from the time variable t to temperature or, equiva- 
lently, the variable x. 

For x above the freezing point xf, Y » Y eq and, upon omitting the contributions of the derivative terms dh/dx, an 
approximation which is also adopted in the standard cosmological treatments ^3 > we obtain for the solution of l|3.4|l 



Y eq = exp 



TH- 1 



-dx 



Here, Y^P corresponds to n£q and in the non-relativistic limit is given by 



Y. 



(0) = ^-3/2, 

eq 2ir 2 h 



21 ( 2TO )-^ exp 



(3.10) 



(3.11) 



where g s counts the particle's spin degrees of freedom. 

In the regime x < Xf, Y >> Y^ the equation (|3.4|l can be written as 

^7 = ~ m x ( w ) ( —G N geff ) h + —y- (3.12) 

Applying (|3.12() at the freezing point Xf and using (|3.10() and (|3.11|l . leads, after a straightforward calculation, to the 
determination of Xf = Tf / through 



Xf 1 = In 



0.03824 g t 



Mpi anc k m<i 1/2 



1 , 9* 

7: in [ — 

2 



rif- 



(3.13) 



As usual, all quantities are expressed in terms of the dimensionless x = T/m^ and Xi n corresponds to the time to 
discussed previously, taken to represent the exit from the inflationary period of the Universe. 

The first term on the right-hand-side of (|3.13|) is that of a conventional Cosmology for, say, an LSP carrying g s spin 
degrees of freedom, playing the role of the dominant Cold Dark Matter species in a concrete and physically promising 
example which we use in this work. The quantity (va) j is the thermal average of va at Xf and is g e ff of 
conventional Cosmology at the freeze-out point. The same notation holds for g*. In our treatment above, we chose 
in H3.13fl to present Xf in such a way so as to separate the conventional contributions, which reside in the first term, 
from the contributions of the dilaton and the non-critical-string dynamics, which are contained within the last two 
terms. The latter induce a shift in the freeze-out temperature. The penultimate term on the right hand side of l|3.13|l . 
due to its logarithmic nature, does not affect much the freeze-out temperature. The last term, on the other hand, is 
more important and, depending on its sign, may shift the freeze-out point to earlier or later times. To quantify the 
amount of the shift one must solve the equations l|2.11|l of ref . |l| • We shall do this in a forthcoming publication , 
where we shall also present a detailed analysis of the effects of the non-critical string and dilaton-source terms on the 
constraints on supersymmetric particle physics models, extending conventional cosmology works Q . 

For our purposes here we note that, in order to calculate the relic abundance, we must solve (I3.12|) from Xf to 
today's value xq, corresponding to a temperature Tq ~ 2.7°K. Following the usual approximations we arrive at the 
result: 



Y-\x )=Y- 



l (xf) + (^r) TUX M Planck g* 



h(x )J 



Xf TH -1 



xY 



-dx 



(3.14) 



In conventional Cosmology ^3] 3* is replaced by g+ and the last term in Ij3.14|l is absent. The quantity J is J = 
Sxo ( v<7 ^ ^ replacing Y(xf) by its equilibrium value l|3.10|l the ratio of the first term on the r.h.s. of (|3.14|l to 
the second is found to be exactly the same as in the no-dilaton case. Therefore, by the same token as in conventional 
Cosmology, the first term can be safely omitted, as long as 1/ is of order of 1/10 or less. Furthermore, the integral on 
the r.h.s. of l|3.14|l can be simplified if one uses the fact that (va) n is small as compared with the expansion rate a/ a 
after decoupling . For the purposes of the evaluation of this integral, therefore, this term can be omitted in (|3.12|) . as 
long as we stay within the decoupling regime, and one obtains: 



d_l_ 
~dx~Y 



TH- 



xY 



(3.15) 
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By integration this yields Y(x) = Y(xq) exp( — J TH 1 dx/x ). Using this inside the integral in l|3.14[l we get 



I / ft TH- 1 , \ 1 , 7T 



(h(xo)Y(x )) = [ 1 + j ~^) dx ) (45) m * M ^„ cfc 5* 5 J (3.16) 

where the function ^p(x) is given by tp(x) = a;exp( — f TH~ 1 dx/x ). With the exception of the prefactor on the 
r.h.s. of 13.16f) . this is identical in form to the result derived in standard treatments, if is replaced by g+ and the 
value of Xf, implicitly involved in the integral J, is replaced by its value found in ordinary treatments in which the 
dilaton-dynamics and non-critical-string effects are absent. 
The matter density of species x is then given by 

Px 1 \ 45 J \tJ M Planck J 1 ' 



where the prefactor / is: 



/ = 1 



x f TH- 1 

x 1p(x) 



It is important to recall that the thermal degrees of freedom are counted by g e // (c.f. H3.6fl ). and not <7e//, the latter 
being merely a convenient device connecting the total energy, thermal and non-thermal, to the temperature T (c.f. 
(|3~5|) ~). Hence, 

Tj\ 3 _ 9eff (lMeV) 4 _ 43 1 (3 ig) 



T i) 9eff( T x) 11 11 9* 

In deriving l|3.18|l only the thermal content of the Universe is used, while the dilaton and the non-critical terms do 
not participate. Therefore the x's matter density is given by 

' X V 45 7 11 M Planck g,J- V ' ' 

This formula tacitly assumes that the x s decoupled before neutrinos. For the relic abundance, then, we derive the 
following approximate result 

Ogfcg = (%^)„ rc x ( g A 1>2 (1 + n^rdx) . (3.20) 



,9* J V Jx i>{x) 

The quantity referred to as no-source is the well known no-source expression 

= L ° 66 X GCY " (3.21) 

V X 0) n o-source M P i anck ^ J 

where J = f s (vo~) dx. However, as already remarked, the end point Xf in the integration is the shifted freeze-out 
point as determined by Eq. H3.13|l . The merit of casting the relic density in such a form is that it clearly exhibits the 
effect of the presence of the source. Certainly, if an accurate result is required, one can proceed without approximations 
and handle the problem numerically as in the standard treatments. We shall present such a more complete analysis 
in a forthcoming publication |Tfij . 



IV. DISCUSSION: SOURCE EFFECTS ON PARTICLE PHENOMENOLOGY 

From the expression (|3.20(l above, it becomes evident that the effects of the integral involving the source T on the 
relic density evolution may be quite important. Indeed, if T is kept negative at all times, this results in reduction of 
the relic density with time, contrary to what happens in the case where T is positive. In the former case, predictions 
for supersymmetric models |8j can be drastically altered, since the parameter space is enlarged, leaving more room for 
supersymmetry, probably beyond the reach of LHC, even in the case of constrained minimal supersymmetric standard 
models with compact parameter spaces of the embedding minimal supergravity theory . The opposite happens in the 
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case of positive T, where the parameter space is shrunk and predictions can be very restrictive to almost excluding 
supersymmetry, especially if the prefactor turns out to be a large number. 

In order to get a rough picture of the importance of such changes in the calculation of the relic density, let us 
assume that a/ a + $ = 0, as would be the case in the present era if matter were absent jH]. The non-critical terms 
would also contribute little in this case and from i|3.1[l one would have T w $. Therefore, TH" 1 —1 and the 
function ip{ x ) would be ip(x) x 2 /xq- As a result, the prefactor in H3.2()(l becomes xo/xf which is an enormously 
small number, of order ~ 10 -13 or so. Such a small number would result practically to no cosmological constraints 
on supersymmctric models ! However this situation is not realized in nature, since at the present era neither matter 
nor the non-critical term contributions are negligible 0. A decent approach is to solve the cosmological equations 
(|2.11|) and thus obtain the function T(t) at all times, t, from today to the remote past. This would allow for precise 
information to be obtained on the value of the freezing point Xf and of the prefactor appearing in the relic density 
(|3.2U|) . Thus, a complete numerical treatment, along the lines presented in 0, is needed in order to tackle important 
phcnomcnological questions regarding the relic abundance of Dark Matter in this framework. 

Nevertheless, the smallness of the prefactor in the crude approximation we have employed above, gives us a sneak 
preview of the drastic changes one may be faced with, as a result of the non-critical string dynamics. Moreover, in such 
a framework one could also tackle other important issues, such as a possible resolution of the gravitino overproduction 
problem in effective supergravity inflationary models In the non-critical string framework discussed above, the 
gravitino, as a member of the gravitational supermultiplet of the string, will feel the non-critical and dilaton dissipation 
terms, while matter effects are unimportant for this issue. Thus, its density may be substantially reduced, in a way 
similar to that discussed above, as compared with conventional supergravity scenarios [TtI . affair a suivre. . . 
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